Multiwavelength Observations of a Flare from Markarian 501 
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ABSTRACT 



We present multiwavelength observations of the BL Lacertae object Markarian 501 
(Mrk 501) in 1997 between April 8 and April 19. Evidence of correlated variability 
is seen in very high energy (VHE, E ^ 350 GeV) 7-ray observations taken with 
the Whipple Observatory 7-ray telescope, data from the Oriented Scintillation 
Spectrometer Experiment (OSSE) of the Compton Gamma-Ray Observatory (CGRO), 

■ and quicklook results from the All-Sky Monitor (ASM) of the Rossi X-ray Timing 
Explorer (RXTE) while the Energetic Gamma-Ray Experiment Telescope (EGRET) 

. did not detect Mrk 501. Short term optical correlations are not conclusive but the 

U-band flux observed with the 1.2m telescope of the Whipple Observatory was 10% 
\ higher than in March. The average energy output of Mrk 501 appears to peak in the 2 

"*^h \ keV to 100 keV range suggesting an extension of the synchrotron emission to at least 

Qlh| 100 keV, the highest observed in a blazar and ~100 times higher than that seen in the 

other TeV-emitting BL Lac object, Mrk 421. The VHE 7-ray flux observed during this 

■ period is the highest ever detected from this object. The VHE 7-ray energy output is 
somewhat lower than the 2-100 keV range but the variability amplitude is larger. The 

£> . correlations seen here do not require relativistic beaming of the emission unless the 

^ , VHE spectrum extends to ;>5 TeV. 
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Subject headings: BL Lacertae objects: individual (Markarian 501) — gamma rays: 
observations 



1. Introduction 



Markarian 501 (Mrk 501) is a BL Lacertae object (BL Lac) at z = 0.033 (e.g., Stickel, Fried 
fc Kiihr 1993j ) making it the second closest known BL Lac after Mrk 421. Mrk 501 was discovered 



as a 7-ray source at E > 300 GeV by the Whipple collaboration ( Quinn et al. 1996 ) and recently 



confirmed by the HEGRA collaboration ( Bradbury et al. 1997 ). Its initial detection level was 8% 
of the VHE flux from the Crab Nebula. Like the other TeV-emitting AGN, Mrk 421, Mrk 501 
exhibits day-scale variability in its VHE 7-ray emission ( |Quinn et al. 19*96 ), but prior to 1997, 



the 7-ray emission from Mrk 421 was generally observed to have a higher mean flux (Schubnell et 
al. 1996| ), to be more frequently variable ( |Buckley et al. 1996] ), and exhibit higher amplitude and 
shorter time-scale flares (Gaidos et al. 1996). EGRET has never detected Mrk 501 ( |Fichtel et al 



1994] |Thompson et al. 1995| ). Its low VHE flux and the nearby EGRET source, 4C 38.41 ( Mattox 



et al. 1993| ), were previously assumed to account for this. 



Here, we present Whipple Observatory VHE 7-ray and optical, OSSE, and EGRET 
observations of Mrk 501 during a flare in 1997 April. We also include quicklook X-ray results from 
the ASM. Overlapping observations of Mrk 501 with RXTE QCatanese 1997|) and Beppo SAX 
( Pian et al. 1997]) are reported elsewhere. Details of the observations are given in § 2. They show 
evidence of a correlation among the wavebands covered by the above experiments (see § 3), except 
EGRET which detected no emission from Mrk 501. The implications of these results are discussed 
4. 
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2. Observations and Analysis 

2.1. Whipple Observatory observations 

The VHE observations reported in this paper were made with the atmospheric Gerenkov 
imaging technique (Cawley &: Weekes 1995) using the 10 m optical reflector located at the 



Whipple Observatory on Mt. Hopkins in Arizona (elevation 2.3 km) ( pawley et al. 1990| ). A 
151 photomultiplier tube (upgraded in 1996 from 109) camera mounted in the focal plane of the 
reflector records images of atmospheric Cerenkov radiation from air showers produced by 7-rays 
and cosmic rays. The 7-ray selection was based on the Supercuts criteria (Reynolds et al. 1993), 
however the camera upgrade necessitated that the analysis cuts be changed slightly ( |Weekes et al 
1997 ). The data reported in this paper are analyzed with a Tracking analysis (cf. |Quinn et al 



1996). 
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Count rates for 7-rays are converted to integral fluxes by expressing them as a multiple 
of the Crab Nebula count rate and then multiplying that fraction by the Crab Nebula flux: 
/(> 350 GeV) = (1.05 ± 0.24) x 10~ 10 photons cm" 2 s" 1 (|Hillas et al. 1997| ). The errors in the 



absolute fluxes listed below are in general dominated by the uncertainty in the Crab Nebula flux. 
The statistical significance of the detections are more accurately indicated by the 7-ray count rates 
(see Fig. |l]a). This procedure assumes that the VHE 7-ray flux of the Crab Nebula is steady, as 7 
years of Whipple Observatory data indicate ( Hillas et al. 1997| ), and that the Mrk 501 spectrum is 



identical to that of the Crab Nebula between 0.3 and 10 TeV, dN/dE oc E~ 2 - 5 ( [Hillas et al. 1997[ ) 



which may not be the case. Estimates of the VHE spectrum of Markarian 501 will be undertaken 
in future work. 

VHE observations of Mrk 501 were taken nightly April 7 - 19, for a total exposure of 19.9 
hours. The mean flux above 350 GeV during this period was (16.4±3.9)xl0~ n photons cm s , 
which is ~1.6 times the flux of the Crab Nebula. By comparison, the 7-ray flux from Mrk 501 was 



never recorded to be higher than 0.5 times the Crab Nebula flux in previous seasons (Quinn et al 
1996). The flux during April ranged from a low of (4.9±1.8) x 10 -11 photons cm~ 2 s -1 on April 
19 to a high of (40.5±9.6) x 10~ n photons cm 2 s 1 on April 16. The latter is the highest flux 
recorded from Mrk 501 by the Whipple Observatory. No evidence of hour-scale variability was 
found within this data set. 

Optical observations were taken April 7-15 with the 1.2m Ritchey-Chretien telescope of the 
Whipple Observatory using standard UBVRI filters. The 12' full field CCD frames were analyzed 
using relative photometry with a 6" aperture applied to Mrk 501 and several comparison stars in 
the field of view. The U-band fluxes (see Fig. |l|d) are expressed in arbitrary units with no galaxy 
light subtraction. These results are described in detail by Buckley & McEnery (1997). 



2.2. Compton Gamma-Ray Observatory Observations 

CGRO observed Mrk 501 for the period April 9 - 15 as the result of a Target of Opportunity 



initiated in response to the reported high level of VHE 7-ray activity (Breslin et al. 1997). We 



report here on the measurements of the OSSE and EGRET instruments on CGRO. 

OSSE observes 7-rays in the 0.05 - 10 MeV energy range. A description of OSSE, its 
performance and analysis procedures is given by Johnson et al. (1993). For the observation 
of Mrk 501, the CGRO viewing attitude and the OSSE observing sequence were selected to 
avoid the X-ray binary Hercules X-l and the EGRET source 4C 38.41. This was achieved by 
viewing Mrk 501 ~ 0?6 off axis, or at ~ 80% of maximum response, and by making background 
observations which were offset by ±7?5. 

The OSSE data taken during this 6 day period yielded a very strong detection of Mrk 501 
(see Table ||). The daily flux observed by OSSE varied by over a factor of 2, peaking on April 13 
(see Fig. ||). Mrk 501 has not previously been observed by OSSE, so no comparison to previous 
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flux levels can be made. The average 50 — 150 keV flux detected from Mrk 501 is higher than the 
emission seen by OSSE from any other blazar ( [McNaron-Brown et al. 1995 ) except the highest 



recorded state from 3C 273 ( [McNaron-Brown et al. 1996 ). We note that it may not be strictly 



valid to compare Mrk 501 to 3C 273 because the latter object is not a typical blazar. The average 
of the OSSE data for this 6 day observation is well fit by a single power law in the energy range 
0.05 - 2.0 MeV (x 2 /dof= 0-794, 64 dof) of the form: 

/ E \ - 2 - 08±0 - 15 
N(£7) = (1.8 ± 0.1) x 10~ 2 — — — ) photons cm" 2 s _1 MeV -1 



O.lMeV; 

where the uncertainties represent the 68% confidence interval for joint variation of two parameters 
(Ax 2 = 2.3). The spectral index is somewhat steep compared to OSSE spectra of other blazars 
but it is still within the 1.0 - 2.1 range reported by McNaron-Brown et al. (1995). There was no 
statistically significant variation in the power law index with time when fitting power law spectra 
to the 6 daily average OSSE measurements. 



The EGRET instrument (e.g., Kanbach et al. 1988; Thompson et al. 1993j ) is sensitive to 



7-rays in the energy range spanning approximately 30 MeV - 30 GeV. The EGRET data were 



analyzed with maximum likelihood techniques ( Mattox et al. 1996 ), taking into account the 



contributions of isotropic and Galactic diffuse emission ( [Hunter et al. 1997 ) and other bright 



sources near the position of Mrk 501. The observations taken April 9-15 indicate an excess of 
1.5cr, which corresponds to a 2a upper limit of /(>100 MeV) < 3.6 xlO photons cm~ 2 s . 



3. Results 

Figure |l] shows daily flux levels for the contemporaneous observations of Mrk 501. The 
average flux level in the U-band in March is also included in the figure (dashed line in Fig. |]d) to 
indicate the significant (^10%) increase in flux between March and April. An 11 day rise and fall 
in flux is evident in the VHE and ASM wavebands with peaks on April 13 and 16. The 50-150 keV 
flux detected by OSSE also increases between April 9 and 15 with a peak on April 13. The optical 
data may show a correlated rise but the variation is small (at most 6%). Subtraction of the galaxy 
light contribution will increase the amplitude of this variation, but it should still remain lower 
than in X-rays given that the R-band contribution of the galaxy light is ~75% (Wurtz, Stocke, 



fc Yee 1996| ) and the U-band contribution should be much less. The ratio of the fluxes between 



April 13 and April 9 are 4.2, 2.6, 1.7, and 1.01 for the VHE, OSSE, ASM, and U-band emission 
respectively. 

In Figure || we show the spectral energy distribution, plotted as the power per logarithmic 
bandwidth, vF u , versus frequency, v. The average flux for April 9 - 15, the duration of the CGRO 
ToO, is indicated by the filled circles. The filled square indicates the peak VHE 7-ray flux from 
April 16. The OSSE data for this plot have been divided into 8 energy bands between 50 keV and 
470 keV. No excess emission is seen above 470 keV so those data are not shown. The VHE integral 
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fluxes are converted to single point fluxes by assuming the spectrum follows a power law with the 
same spectral index as the Crab Nebula. The EGRET integral flux upper limit is converted to a 
single point flux limit assuming dN/dE oc E~ 2 . The ASM 2-10 keV count rate is converted to a 
flux point by normalizing to the Crab flux and assuming Mrk 501 has the same photon spectrum 
as the Crab (E~ 2m ) within that energy range. From Fig. a peak in the power output occurs 
within the 2-100 keV range and the VHE 7-ray power output is much higher than at EGRET 
energies but somewhat below that in the 2-100 keV range. 



4. Discussion 



The results presented here show that for Mrk 501, like Mrk 421 ( Buckley et al. 1996 ; 



Macomb et al. 1995 ), the VHE 7-rays and the soft X-rays vary together, the energy budget at 
X-rays and 7-rays is comparable, there is evidence of correlated optical variability, and the cut-off 
in the synchrotron emission occurs at much higher energies than in radio-selected BL Lacs or 
flat-spectrum radio quasars. However, Mrk 421 has a large deficit in energy output between 1 keV 
and 50 keV ( [McNaron-Brown et al. 1995| ; [Buckley et al. 19961 ). With Mrk 501, there is instead a 



peak in the power output at ~100 keV, and the deficit has shifted to between 1 MeV and 350 GeV. 
Also, the variability amplitude for the VHE 7-rays is larger than the X-ray and OSSE variations, 
whereas similar flares in Mrk 421 showed comparable amplitude variations ( Buckley et al. 1996] ; 



Macomb et al. 1995). 



A likely explanation of the OSSE detection is that the synchrotron emission in Mrk 501 
extends to 100 keV, compared with the ~1 keV cutoff seen in Mrk 421. If this is synchrotron 
emission from nonthermal electrons, then the maximum electron Lorentz factor is given by 
7 max = 3 x 10 6 [£' syn (100 keV)]/^] 1 / 2 , where B is the magnetic field strength in Gauss, 8 is the 
Doppler factor of the radiating plasma, and £ , syn (100 keV) is the cutoff energy of the synchrotron 
spectrum in units of 100 keV. The correlated variability between the hard X-rays and the VHE 
7-rays implies that the same nonthermal electrons are responsible for the two emission components 
(cf. IMacomb et al. 1995Q , with the VHE 7-rays produced by Compton scattering of either internal 



synchrotron photons (e.g., [Konigl 1981 ; Maraschi, Ghisellini, & Celotti 1992; Bloom & Marscher 



1996) or photons produced external to the jet (e.g., Dermer, Schickeiser, & Mastichiadis 1992; 



Blandford 1993j ; Sikora, Begelman, & Rees 1994). For low redshift objects, electrons must have 



Lorentz factors 7 > 2 x 10 6 -E'c(TeV)/5 to Compton-scatter photons to energies of -E'c(TeV) TeV. 
Thus the magnetic field B £ 2£ syn (100 keV)<5/[£c(TeV)] 2 (see also |Buckley et al. 1997|) . 



A lower limit to the magnetic field is implied by the variability of the synchrotron 
emission. Electrons cool through synchrotron losses in the co-moving frame on a time scale 
At syn = 7.7 x 10 8 /(£ 2 7) s . The source might intrinsically vary on a longer time scale, so a lower 
limit to the observed variability time-scale for low-redshift objects is given by Ai b s ^ At syn /5. 
This implies that B £ [9000/ (-f5At ohs (days)] 1 / 2 . This expression can be written in terms of the 
observed photon energy i? b s (eV) using the relation v \> s = 2.8 x 10 6 -E>7 2 <5 from synchrotron theory. 
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Putting the above constraints together, the allowed magnetic field strength of the emission region 
is given by 

[ J E obs (eV)<5A4 s (days)]- 1 / 3 £ B & 2£ syn (100 keV)5/[^ c (TeV)] 2 . (1) 

The ASM (£obs = 2000 eV) flux varies on time-scales at least as short as one day, implying that 
B £ 0.086-V 3 G. Such a range of field strengths are expected at the base of a jet from independent 
considerations (e.g., Marschcr fc Gear 1985| ). 



Equation |T| provides a new test for beaming in 7-ray blazars, although it assumes that the 
high-energy component arises from Compton-scattering of photons by the same nonthermal 
electrons which produce the synchrotron component. If the lower limit in equation |l| exceeds the 
upper limit when S is set equal to 1, then beaming is implied. Here it is especially crucial to 
determine the cutoff energy at TeV energies. If the high-energy spectrum of Mrk 501 continues 
without break to ^ 5 TeV, then beaming is implied. 

This beaming test complements 77 transparency arguments which also provide lower limits 
to 5. The > 350 GeV 7-rays interact with the X-ray photons in the extreme relativistic regime 
of pair production (i.e., -EVhe-E-x m^c 4 ), which has been treated by Dermer & Gehrels (1995). 
One finds that t 77 <^ 0.1 for this interaction, so no beaming is required. Improved limits to the 
Doppler factor through 77 transparency arguments and from equation [l] will be obtained if short 
time scale correlations are detected between VHE 7-rays and photons in the EUVE or optical 
regime. Also, if the spectrum observed with OSSE extends beyond 0.511 MeV, 77 transparency 
arguments can be applied to this 7-ray emission directly ( McNaron-Brown et al. 1995| ). 



Finally, it will be important to investigate whether the extension of the synchrotron emission 
in Mrk 501 to 100 keV is always present, or only occurred during the flaring activity shown in this 
work. Previous compilations of X-ray spectra for Mrk 501 spanning 0.03 - 30 keV (e.g., |Ciliegi et 



al. 1993j ) are not conclusive, mainly due to the low statistics at the high energies. Further studies 



of Mrk 501 with OSSE and the high energy instruments on RXTE and BeppoSAX could resolve 
this issue. 

We acknowledge the technical assistance of K. Harris. This research is supported by grants 
from the U. S. Department of Energy and by NASA, by PPARC in the UK and by Forbairt in 
Ireland. The X-ray data in this work are quicklook results presented by the ASM/RXTE team. 
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Table 1. OSSE flux measurements for Mrk 501 



0.05-0.15 MeV 


0.15-0.50 MeV 


0.5-3.5 MeV 


Flux a 24.2±1.3 


2.0±0.3 


<0.2 



a Flux units are 10 3 photons cm 2 s 1 MeV l . The 
upper limit is 2a. 
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Fig. 1.— (a) VHE 7-ray, (b) OSSE 50 keV - 150 keV, (c) ASM 2-10 keV, and (d) U-band optical 
light-curves of Mrk 501 for the period 1997 April 2 (MJD 50540) to 1997 April 20 (MJD 50558). 
The dashed line in (d) indicates the average U-band flux in 1997 March. 
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Fig. 2. — The spectral energy distribution of Mrk 501. The contemporaneous observations taken 
April 9-15 (filled circles), archival data (open circles) ( |Gear et al. 1994 ; pwen et al. 198C ; Impey 
fc Neugebauer 1988| ; [Kotilainen et al. 1992fc |Edelson et al. 199$ [Singh fe Garmirej ; |Staubert et 



al. 1986| ; |Fichtel et al. 1994| ; |Bradbury et al. 1997| ), the mean VHE 7-ray flux in 1995 (filled star) 
( Quinn et al. 1996| ), the mean flux in 1996 (open star), and the maximum VHE 7-ray flux, detected 
on April 16, (filled square) are indicated. 



